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Using the dynamic renormalization group (DRG) tech- 
nique, we analyze general nonlinearities in a conservative non- 
linear growth equation with non-conserved gaussian white 
noise. We show that they fall in two classes only: the 
Edwards- Wilkinson and Lai-Das Sarma types, by explicitly 
computing the associated amputated two and three point 
functions at the first order in perturbation parameter(s). We 
further generalize this analysis to higher order nonlinearities 
and also suggest a physically meaningful geometric interpre- 
tation of the same. 



Growth of thin films from vapor has been a subject 
of significant interest for both experimentalists and the- 
orists. The non-equilibrium nature of this growth makes 
the techniques applicable to the study of fractals, avail- 
able for understanding the physics of thin film growth; 
particularly in the thermodynamic limit |jj . Most of the 
characterization can be done by observing spatio tem- 
poral behavior of the interface fluctuations. This be- 
havior can be studied within the framework of contin- 
uum equations such as a Langevin type equation pro- 
posed by Edwards and Wilkinson (EW) 0. It allows a 
non-equilibrium statistical mechanical description which 
brings out the scale invariant nature of the interface 
roughness. Consequently one obtains two dynamical crit- 
ical exponents a and [3 (or equivalently a and z — a/f3) 
which completely characterize the growth. The corre- 
lation length along the substrate varies as r a and t x / z 
whereas that along the normal to the surface exhibits a 

behavior. 

The geometrical feature of growth represented by the 
lateral component of the growth velocity was captured 
by Kardar, Parisi and Zhang (KPZ) ||, by introducing 
a nonlinear term proportional to (Vh) 2 in the EW equa- 
tion: 



dh 
dt 



v\7 2 h+^(Vhf 



V 



(1) 



where h(x,t) is the height function and r)(x,t) is a non- 
conserved gaussian white noise characterized completely 
by its one and two-point functions: 

fa(x,i)>=0 

and (r?(x, t)r)(x',t')) = 2D5(x - x')S(t - t') (2) 



Dynamic renormalization group technique Jfj applied 
to the KPZ model yields z = 3/2 and a — 1/2, based 
on one-loop calculations for a one dimensional substrate. 
These values of a and z are obtained by looking for a non- 
trivial fixed point for the set of one-loop flow equations 
for the parameters v, A and D stated above. It turns out 
that these results are exact in the sense that they hold to 
all orders and can be obtained by invoking the Galilean 
invariance of the KPZ equation and the Fokker-Planck 
equation associated with it for the 1 + 1 case jlj. The 
KPZ equation, however, is of the non-conservative type, 
(i.e. it can not be written in the form ^ = V • j where 
j(x, t) is a local current) and so is relevant mainly for 
growth with overhangs and vacancies corresponding to 
the regime of high temperature and high volatility of the 
deposited material. 

A general fourth order conservative equation was pro- 
posed by Lai and Das Sarma (LDS) || of the form 
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If X 22 = then in the absence of noise, this equation 
is symmetric under h — > — h. When v 2 is zero the V ■ 
(V/1) 3 term is more relevant than the V 2 (V/i) 2 term as 
a perturbation to the linear term V h. If one makes the 
coefficients in this equation time dependent, the situation 
can change significantly with the possibility of the latter 
term taking over the former one. Interpretation of this 
in terms of different relaxation mechanisms operating at 
different times during the pre-asymptotic region in time 
will be discussed elsewhere M. 

In a recent paper, Das Sarma and Kotlyar Q, have 
shown that the most general conserved fourth order non- 
linearity, viz. A13V • (V/1) 3 belongs to the universality 
class of the EW model., i.e. the interaction produces a 
—v 2 k 2 term in the inverse propagator at zero frequency, 
where the coefficient v 2 is proportional to the first power 
of the perturbation parameter A13. This confirmed the 
earlier numerical evidence based on simulations in (1 + 1) 
dimensions due to Kim and Das Sarma 

The purpose of this letter is to show that this analy- 
sis is in fact completely general and a general nonlinear 
conservative growth equation reduces to either Lai-Das 
Sarma or the EW form. We shall consider the following 
equation with a general nonlinearity of order(2n + 2) for 
integer n, 
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Note that the linear term too is of order (2n + 2), in order 
to keep the nonlinear terms relevant in RG sense. 

Consider e.g., the general nonlinearity: V • (Vh) 2n+1 . 
(The choices n = and n = 1 correspond respectively 
to the EW equation and reference [6]). For the purpose 
of diagrammatic perturbation theory, it corresponds to 
a "vertex" with (2n + 1) prongs ending on a cross and 
one free ended prong. (see, e.g., the book by Barabasi 
and Stanley for notations and conventions). The ver- 
tex and the corresponding Feynman rule are in Fig.l. 
for example n = 1 it corresponds to the vertex factor of 
|Ai(qi • q2)k • (k — qi — q2), which upon suitable iden- 
tification of the parameters and a change of variables 
becomes the expression of reference [6]. The propagator 
modification at order Ai is shown in Fig. 2. 

Again, with n = 1 we recover the one-loop Hartree 
term of reference [6] . The n-loop Hartree diagram shown 
in figure 2, is equally easy to evaluate; thanks to the 
peculiar nature of the noise contraction which generates 
the Feynman factor 2D<5(k 1 + k 2 )5(u;i + w 2 ). The inverse 
free propagator is given by 



G \k,uj) = -iuj+{-l) n vk 



2 n+2 



(5) 



For our purpose (which is to identify the modification 
of the equation of motion or the hamiltonian whose h- 
variation produces the equation), it is enough to compute 
the 'vertex function' or the 'amputated' correlation func- 
tion which amounts to snapping off the external propa- 
gators. There is a symmetry factor associated with this 
diagram given by 5 (2n + 1) ■ 2n(2n — 1). Also every noise 
contraction gives a 2D factor yielding an overall mul- 
tiplier as (-1)™+ 1 A 1 • n(2n - l)(2D) n . Using the stan- 
dard techniques to evaluate the frequency-momentum in- 
tegrals, we finally arrive at the amputated n-loop Hartree 
vertex expression as 



(-l)"+ 1 A 1 (£>/i/)"n(2n-l)^ 



d-(2n+l) 



dq 
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where Kd is related to the usual angular integration fac- 
tor as Kd = *\<i Sd- It is clear from this expression that 

the V- (\7h) 2n+1 nonlinearity produces a V 2 /i term (iden- 
tified as the coefficient of — k 2 in the above expression) 
at the same order in Ai in dimensions d < d cr i t . This 
puts the "odd order" nonlinearities in the EW class. 

Exactly similar analysis (see Fig. 3 ) of the "even or- 
der" nonlinearity t^A 2 V 2 (V/i) 2 ™ produces the amputated 
three point function (which is just the object that appears 
in the equation) given by 



fc 2 [q(k-q) 



3A 2 
n 



"(2n+l) dp 
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Again, we simply identify from this expression that the 
Lai-Das Sarma term, V 2 (Vft.) 2 is produced at the first 
order in A 2 in dimensions less than the critical dimen- 
sion, putting the even order nonlinearities in the Lai-Das 
Sarma class. 

These arguments are trivially extendible to the non- 
linearities of the form V 2m+1 • (Vh) 2n+1 and \7 2m (Vh) 2n 
(with the appropriate linear terms in the growth equation 
to keep the nonlinearities relevant). They respectively 
generate the generalized EW and the Lai-Das Sarma 
terms, viz. V 2m V 2 /i and V 2m (Vft,) 2 types. Gener- 
ally a nonlinearity of the kind V 2 [V 2m ft-] 2 ™ collapses to 
V 2 (V 2m /i) 2 and V • [V 2m + 1 /i] 2 "+ 1 to V 2 [V 2m h}. 

The above mechanism of generation of the lower order 
terms via higher order nonlinearities reflects onto a sim- 
ple physically meaningful geometric picture of the evolv- 
ing growth front. Consider e.g., the equation 



dh 



^ = vV 2n + 2 h + \ 1 V-(\7h) 2n+1 +r ] 
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The rate |£ depends upon the variation of (V/i) 2n+1 , 
including its sign and h varies locally to attempt to make 
this rate vanish for a stable growth mode. Note that the 
nature of the response for terms V-(V/i) and V-(V/i) 2n+1 
is exactly the same except for its magnitude, given any 
geometrical configuration on the surface for any value of 
n (i.e. the power of the gradient). This fact seems, ac- 
cording to us, to be the root cause for the generation of 
the EW term since the surface morphology in these cases 
must be comparable, considering the compatibility of the 
responses. 

Similar reasoning can be furnished in the case of 
V 2 (Vh) 2n leading to the V 2 (V/i) 2 term. Here the only 
difference is that the growth rate does not depend on the 
sign of h, thus making the LDS term as the lowest or- 
der term allowed. This leads to the conclusion that in 
general the nonlinearities are characterized by the order 
of the variation of h locally, i.e. slope, curvature, etc., 
corresponding to Vh, V 2 h, ■ ■ ■ etc.. and thus a general 
nonlinearity is of the form V 2m {V n h) 2 . 
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FIG. 1. The vertex for V ■ (V/i) 2 " +1 nonlinearity. The 
associated Feynman rule is: 

( (2n+l) Al (qi ' • • • (q2n-l ■ q2n) [k ■ q 2 „+l](5(k - 

i=i qO^-E^i wo 

FIG. 2. The n-Loop Hartree diagram generating the EW 
term with a coefficient proportional to Ai. Also shown is the 
two-point vertex function. 

FIG. 3. The diagram for the nonlinearity of the type 
V 2 (Vfo) 2 proportional to A2 generated by the nonlinearity 
2^A2V 2 (V/i) 2n . Also shown is the three-point vertex func- 
tion. 
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